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Abstract 



It is proven that the inverse localization length of an Anderson model on a strip of 
width L is bounded above by L/X 2 for small values of the coupling constant A of the 
disordered potential. For this purpose, a formalism is developed in order to calculate the 
bottom Lyapunov exponent associated with random products of large symplectic matrices 
perturbatively in the coupling constant of the randomness. 



1 Main result and discussion 

The Anderson model describes electronic waves scattered at random obstacles. Here the physical 
space is supposed to be quasi one dimensional and given by an infinite strip of finite but 
large width L. Hence the Hilbert space is £ 2 (Z, C L ) and states therein will be decomposed as 
ip = (i[}(n)) n ^z with ip(n) 6 C L . The Anderson Hamiltonian on a strip is then denned by 



Here Ax : C L — > C L is the transverse (one dimensional) discrete laplacian with periodic bound- 
ary conditions; denoting the cyclic shift on C L by S, it is given by Ax, = —S — S*. For 
L — 1,2, one may rather set Ai = and A2 = —S = —S*, but our main interest will be in 
the case L > 3 anyway. Furthermore, A G M is the coupling constant of the random potential 
V(ri) : C L — > C L which is a diagonal matrix V(n) = diag(t>(n, 1), . . . ,v(n,L)). All real num- 
bers (v(n,l))nzz,i=i,...,L are independent and identically distributed centered real variables with 
unit variance. Given a fixed energy E e R, it is convenient to rewrite the eigenvalue equation 
Hl(X)iJ} = Eip using the transfer matrices 



(H L (\)i>)(n) 



tp(n + 1) - ip{n - 1) + A L ip(n) + \V(n)ip(n) . 




(1) 



This work concerns the study of the asymptotics of their random products which is char- 
acterized by the Lyapunov exponents. One way |BL| ICLj to define these exponents is to use a 
formalism of second quantization. For p = 1, . . . , L, let A P C 2L denote the Hilbert space of the 
anti-symmetrized p-fold tensor products of C 2L , the scalar product being given via the deter- 
minant. Given a linear map T on C 2i , its second quantized A P T on A P C 2L is defined as usual. 
Now the whole family of non-negative Lyapunov exponents 71 > 72 > • • • > 7l > are defined 
by: 
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where the expectation is taken over all random variables. Positivity of the bottom Lyapunov 
exponent 7^ was already known |BL| IGMt IGRj and is sufficient to assure Anderson localization 
KLS| IBL] . The object of this work is to provide a quantitative lower bound for small coupling 
of the random potential. 

We will need a hypothesis excluding exceptional energies with Kappus-Wegner breakdown 
of the perturbation theory to leading order |KWj . This can in principle be overcome |KW[ lBK| 
ICKj . Let [c] denote the integer part of c G f . 

Main hypothesis: Let fii = —2 cos(^) — E and, if \fi t \ < 2, set e m = -(//; + z^4 — fj,f). 
Except if a — — 1 and j = k or {k, 1} = {m, j}, we suppose 



^ 1 



^ 1 



k,l,m,j = 1, . . . 



L + l 



a = ±l. (3) 



Theorem 1 Suppose E el is in the spectrum of H L (0) and satisfies the Main hypothesis. 

(i) Then 

JL > ^ + 0(X 3 ) . 

(ii) Let Eb be a band edge of Hl(0) and E = Eb + e be in the spectrum of Hl(0). Then 



For the case L = 1, this was proven by Pastur and Figotin |PF| and in a related situation 
of hamiltonian stochastic differential equations by Arnold, Papanicolaou and Wihstutz |AP Wj . 
Actually, we go beyond the above theorem and prove an asymptotic formula for jl in Theorem 
El in Section E~o1 below. We then argue (non-rigorously) in Section El that the bound (i) gives the 
right order of magnitude for all energies away from the band edges and the so-called internal 
band edges which are defined by the property that rji — for some I. Near the band edges, 
the bottom Lyapunov exponent is much larger according to item (ii). Indications for such 
stronger localization properties in this regime appeared also in Kloj . It is straightforward to 
analyse the large deviations of the growth behavior of the transfer matrices around the typical 



2 



behavior given by the Lyapunov exponent with the techniques of |JSS| Section 5]. If one adds 
the supplementary hypothesis that E(v(n, I) 3 ) = 0, then the corrections are actually of the order 
0{\ i ). The main deficiency of the present work is the lack of control of the error term on the 
strip with L and the energy E. 

The method of proof transposes directly to a more abstract setting of random products of 
symplectic matrices with small coupling, if only one supposes that the modulus one eigenvalues 
of the unperturbed part are non-degenerate (in the language of Section |23J elliptic channels are 
then non-degenerate). In order to deal with the degeneracies appearing in the Anderson model, 
the concrete form of the random perturbation in (0) is however heavily used. We believe that 
hamiltonian stochastic differential equations could also be treated. Preliminary results in this 
framework were obtained by Teichert |Teij . 

Let us briefly describe the key steps of the proof. First the transfer matrix at A = is 
diagonalized into symplectic blocks given by rotations (Sections 12.11 and 12. 5|) . Then the matrix 
elements of the random perturbation are calculated in that representation ( Section I2.2j) . This 
normal form allows to derive a basic perturbative formula for the Lyapunov exponents (Sections 
14.21 and 14.6)1 . A new ingredient herein is the consistent use of symplectic frames. It is then 
possible to apply a crucial identity related to the geometry of Lagrangian manifolds (Lemma 
El in Section l3~Tjl . The normal form of the transfer matrix now allows to efficiently control the 
oscillatory sums appearing in the perturbative formula ( Section 14.3)1 . There is an inessential 
technical difficulty due to the presence of so-called hyperbolic channels. They do not appear if 
L is odd and E is near the band center. The text is written such that the reader can understand 
this case and hence the main point of the argument by skipping Section 13.51 and then omitting 
Sections 14.51 and 14.61 

In order to compare Theorem^ with results in the physics literature, let us interprete L as 
the number of channels of a disordered wire and as the associated inverse localization length. 
Then the behaviour 

1l ~ ^ + 0(X 3 ) (4) 

confirms the predictions of Thouless |Thoj as well as the Dorokhov-Mello-Pereyra-Kumar theory 
(see |Benj for a review on the latter). 

For the Anderson model in two dimensions, all waves are expected to localize even at small 
disorder |AALR| . Few rigorous results indicating such a phenomenon are known. Even though 
one may think of the strip as an approximation to the two dimensional situation, it is unlikely 
that (J1J) gives much insight. For a proof of localization, one would need to prove the so-called 
initial length scale estimate in order to apply the multiscale analysis |FS) IDK) IGK) . It states 
that the wave functions on a square of appropriate diameter decrease from center to boundary 
(or inversely) with a high probability. But even if the error term in Q could be neglected at 
small A, exponential decay of typical eigenfunctions is noticeable only on a length scale N given 
by N'Jl = 0(1), that is N = L/X 2 , which is much larger than the strip width L. Therefore 
(J3J) is of interest only in the quasi one dimensional situation. Indeed, Anderson localization 
in two dimensions is expected to be a non-perturbative phenomenon (like BCS theory) and 



3 



thus not tractable by a "naive perturbation theory" as developed here. This is reflected by the 
prediction that the 2D localization length behaves non-analytically like e 1//A for small A [AALRJ. 
Rigorously known is only a lower bound on the phase-space localization of the eigenfunctions 
|SSW| . 

2 Analysis of the transfer matrix 

Each transfer matrix T(n) is a random element of the symplectic group 



SP(2L,M) = {T eM 2Lx2L (R) \T t JT = j} , 




The aim of this section is to construct a symplectic basis transformation M GSP(2L, R) such 
that 

M~ l T{n)M = R(l-\P(n)), (5) 

where the free transfer matrix R (i.e. A = 0) takes a particularly simple form given by a direct 
sum of rotations. The random perturbation P(n) lies in the Lie algebra sp(2L, R). Its matrix 
elements and some of their expectation values will be calculated below. Throughout this section 
the index n is kept fixed and will thus be suppressed. 



2.1 Normal form without disorder 

Let us introduce, for I = 1, . . . , L, 

( <Mi) \ 



<f>i(k) 



41 



exp 



2irt Ik 



V ML) ) 

Then A^0; = — 2 cos(2tt I / L) 0;. Note that the fundamental = <pi is real; moreover, for even 
L, the vector &l/2 — <Pl/2 is real as well. For other /, real normalized eigenvectors are obtained 
by 



L-lj 



-I 1 



V2 V 1 1 

Next define an orthogonal matrix m G 0(L, R) and unitaries d, f G U(L, C) by 



m = ($i,...,$ L ) 



/ 



m = f d . 



Finally introduce the diagonal matrix fi = diag(/ii, . . . , /i^) where \xi = — 2 cos(27r l/L) — E. 
With these notations, 



m*(A L -£)m = f*(A L — E)f = n . 
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Figure 1: Plot of the the energy levels for L = 13. (i) Here E = 0.95 and = 2. The 
eigenvalues u. , /i 2 , /in and /ii 2 are outside of the window [—2, 2] and hence hyperbolic, (ii) 
Here E = —0.03. All eigenvalues are elliptic. 

An eigenvalue fii will be called elliptic if \fii\ < 2, hyperbolic if > 2 and parabolic 
if = 2. Here energies E for which there are parabolic eigenvalues are excluded by the 
hypothesis (jSJ). An energy E is in the spectrum of Hl(0), the Laplacian on £ 2 (Z, C L ) if there 
exists an elliptic or a parabolic eigenvalue. The spectrum of Hl(0) is hence [—4,4] if L > 2 
is even and [—4,2 + 2 cos(7r/L)] if L > 1 is odd. If L is odd and E slightly above the band 
center (Fig. 1(h)), all eigenvalues are elliptic while for energies E outside of the spectrum, all 
eigenvalues are hyperbolic. In between one has both hyperbolic and elliptic eigenvalues. For 
notational convenience, we will suppose that fjii < 2 for all /. Hence there exists Lh < \ such 
that that fii is hyperbolic for / = 0, . . . , Lh, and elliptic for / = + 1, . . . , W]. Moreover, 
there is a degeneracy Hl-i — due to reflection symmetry which will be further analyzed in 
Section I2~3l In case there are no hyperbolic eigenvalues, let us set = — 1. In Fig. 1 are shown 
examples of (i) a situation with mixed elliptic and hyperbolic eigenvalues and (ii) a situation 
with only elliptic ones. 

For later use, let us set g = diag(gi, . . . , gi) where gi — 1 if /i/ is elliptic and gi — i if \i\ is 
hyperbolic. Note that [d, //] = [d,g] = [fJ,,g] = 0. 

In order to diagonalize the transfer matrix, let us introduce the diagonal complex L x L 
matrix k=| + w^- — las well as the (possibly complex- valued) matrices 



where here and below all roots are taken on the hrst branch. Then one immediately verifies 
that 
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atM Ai r E -})n ' k 



Finally the r.h.s. will be transformed into a normal form which is real-valued and symplectic. 
For elliptic /i/ define < rji < n by k\ = e lVl , for a hyperbolic one t]i > by k\ = e m . Then 
introduce the rotation and hyperbolic rotation matrices by (rj 6 M): 

cos(?7) — sin(?7) \ , . / cosh(?7) sinh(^) 

Setting 



( sin(?7) cos(r]) J ' ^hvJ) I s j n ] 1 (^) cosh(?7) 



one verifies 



G- X C- X ( l !jcG=%)ffi...®%), 

where Ri(t)i) is either R e (r]i) or Rh{rji) depending on whether is elliptic or hyperbolic, and 
where the direct sum is understood such that Ri{rji) acts on the Zth and (I + L)th component 
of C 2L , namely the normal form R = R±(i]i) © ... © Rl{vl) is a real- valued symplectic matrix. 
In case that there are parabolic eigenvalues, the normal form contains Jordan blocs just as in 

Let us resume the above. Set M = NCG, then 

M" 1 ( AL ~ E "q 1 ) M = R. 
Furthermore the symplectic basis change M e SP(2L, M) is more explicitly given by 

/ mh \ / h- l m* 

M = J , M 1 = 

\ m|(fi+-)/i mh' 1 g 2 J \ — \ g 2 (« + ^) /i m* g 2 hm* 

where 



'^T = diag(/i X , ...,h L ) 



K — — 




/i; elliptic, 
/ii hyperbolic. 
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2.2 Calculation of the perturbation 



First note that 



T 



A L -E -1 
1 



1 
1 



A 



Hence P e sp(2L, K) in equation © is given by 
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M 
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g 2 hd*Vdh 



(6) 



Here the identity m = fd was used as well as the definition V = f*Vf. Actually / is the matrix 
of the discrete Fourier transform so that V = V* is the Toeplitz matrix associated with the 
Fourier transform of the random potential (at fixed height n) . More precisely, define 



v(k) 



- «(o ex p ( 
i=i ^ 



2m Ik 
L 



One has v(k) = v(—k) and v(k + L) = v(k). Then 

0(1) 

V 



\v(-L + l) 



v(L-l) \ 

v(l) 
v(0) ) 



2.3 Symplectic channels 

The symplectic channels of a symplectic matrix are by definition the eigenspaces of its (possibly 
degenerate) eigenvalue pairs (k, -). For the elliptic-hyperbolic rotation R = -Ri(^i) © ... © 
Rl{vl), they are the maximal subspaces of C 2i characterized by the property that R rotates 
by the same angle. For the Anderson model on a strip studied here, there are L c + 1 = [^] + 1 
channels. The 0th channel is associated with the fundamental fiL and given by the span of 
{^l,G2l} where ei G C 2L has a non-vanshing entry equal to 1 only in the Ith component. If L 
is even, also the -|th channel is simple and spanned by {ei/2, £31/2}- l^ue to the degeneracy 
Hi = hl-i of the spectrum of A^, all other channels are doubly degenerate. For / 7^ 0, -|, the 
Zth channel is spanned by {ei, ei-i, ei+L, ^2L-i\- Let us denote the degeneracy of the /th channel 
by v\. In accordance with the above, the /th channel is called elliptic if < 2 and hyperbolic 
if \ni I > 2. With these notations, C 2i is hence decomposed into a direct sum of the L c + 1 
channels. 

There are Lh + 1 hyperbolic channels. As we supposed (for convenience) that \i\ < 2, the 
hyperbolic channels are the first ones, namely channels 0, . . . , Lh (compare Fig. l(i)). 
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The projection on the Zth channel will be denoted by 717. It satisfies [tti, J] = 0. The non-zero 
eigenvalues of n\Rn\ are e ±ir)l if the Zth channel is elliptic and e ±m if it is hyperbolic. Hence one 
can decompose 717 into the corresponding eigenspaces 717 = + nf. Resuming all the above, 



i c2L = JJ + n > R < = e±l9im *f ■ (7) 
1=0 

2.4 A privileged basis 

For explicit calculations in the next section, it will be convenient to dispose of a basis of C 2L 
with the following properties: 

1. the basis vectors are orthonormal and compatible with the symplectic structure; 

2. the basis vectors are eigenvectors of R; 

3. the matrix elements of the perturbation P w.r.t. the basis are particularly simple, i.e., up 
to a constant, given by the Fourier transform v of the potential. 

For this purpose let us note that, for any n G K, linearly independent eigenvectors of R e (v) 

are ( \ ) 1 ( ; ), wMle for Z ca„ cll00S e ( \ ) and (\). « 

with equation (jfi]L the third property is verified if the top L components of the basis vectors 
commute with h and are constructed with the inverse d* of d or its quasi- inverse d l . Indeed it 
can be verified that dd l = S is the reflection in C L sending component Z to component (L — l). In 
particular, the components |f and L are left invariant by S. Moreover, S preserves the channels, 

namely [S, d] = [S, g] = [S, h] = and ( ^ ^ j commutes with ni. Hence one is led to define 
the following: 

W = (w+ ...,w+ wr,...,wl) = ( _f gfp ) . (8) 

It can readily be verified that 

WW = 1, W*JW = U r 9 + 9 L i9 *~ 9) ? )■ (9) 

2\(g- g*)S -g-g* J K 1 

Hence the vectors |w^), . . . , |w£), |w^), . . . , |w£) defined by (jSJ) form an orthonormal basis of 
C 2L (we switch to Dirac notation for vectors here). Let us write out examples for a simple and 
a double eigenvalue and fxi, 1 < Z < -|, respectively) more explicitly: 

1 1 

|w£) = -j= (\e L ) -ig L \e 2 L}) , |wf) = - (z |e*} + je^-j) + # |e i+J ) - «^ |e 2 i,-/» . 
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Note that for an elliptic eigenvalue fii, \w t ) is the complex conjugate of |w ; + ), but if \i\ is 
hyperbolic, this only holds for the first L components. As requested, one reads off from Q that 

{± i\wf) iii elliptic , 

(10) 
=F Wi-ii hyperbolic . 

Also the second required property follows from (jHJ: 

R \wf) = e ±imn ' |w±) . 

Finally, 



*o = I w l)( w lU = l w f)( w f I + l w i-;)( w L-ih *t = |w|>(w|| . (11 

2 2 2 



2.5 Matrix elements of the random perturbation 

In this section, we show that the third desired property of the basis stated in Section 12.41 is 
fulfilled and then exploit it in order to calculate the matrix elements of the perturbation and 
some expectation values thereof. Taking into account that [d, h] = 0, it now follows from (jOJ) 
and (H) that 



W*PW - ( 9kVk 9hVhS \ w*P*W - % ( ~ hVh 9* hVhg*S 



2 V SghVh -SghVhS J ' 2 V -ShVhg* ShVhg*S 

Hence one can read off, for two signs r, a, 



(12) 



1 1 

(w[|P|w£) = t -gihih k v(ak-rl) , (w[|P*|w£) = -a -glhih k v(ak - tI) . (13) 
Let us collect some useful identities. 



Lemma 1 Let w ■ G C be unit vectors satisfying 7r • Wj = Wj . Let P = P + P* . 

(i) E (w\P\w'} = for any w,w' G C 2L . 

(ii) If either k ^ I or a ^ r, then 

E (w[\ P \w%) (w%\P\wJ) = -rcr-^- gi g k hfhl . 

(iii) If either k ^ I or o ^ r, then 

E(wJ\P*\wl) (wl\P\wJ) = -Lhftl . 

(iv) Let channels I, k be elliptic. Then rfPiv% = and rfPn^ = 2 TrfPir^ . Moreover, 
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E«|P|tO (w k °\P\w?) = \ h]h\ = E«|P|tO (w?\P\w k °) . 

(v) 

E «| |i?P| 2 K) = ~ /i 2 v fr 2 , W /iere /i 2 v = £ $>* /^ cosh((l - ^) % ) . 

fc=0 

One might think of items (ii) through (v) as follows. Even though the perturbation P lifts 
the degeneracy of the channels, taking expectation values re-establishes it. 

Proof, (i) This follows directly from Et>(/) = for all I. 

(ii) Let wj = aiwj + biw T L _ l with \ai\ 2 + \bi\ 2 = 1 and bi = if the channel is simple, that is 
1 = 0,^. With these notations, it follows from 

(wJ\P\wl){wl\P\wJ) = -oT X - gi g k h]hl- 

■ {alakviak — rl) + ~aibkv(—ak—Tl) + bidkvicrk + rl) + bibkv(—ak+Tl)) ■ 

■ (a/afc-0(— ak+rl) + aibkv{—uk — Tl) + b{akv{ak + Tl) + bibkv(p-k — rl)) . 

Now Ef (Z) 2 = 1 implies that Et>(g){)(j>) = t <5g,-p- With a bit of care, one can now check that 
the expectation value of the product of the last two factors is equal to 1. 

(iii) This is proven in the same manner as (ii). 

(iv) As gi = g k = 1, it follows from (JT2J) that W*PW = W*(P + P*)W has only off-diagonal 
entries in (^ + n k )C 2L equal to twice those of W*PW . The first equality now follows directly 
from (ii), while the second one is checked similarly. 

(v) One first verifies as above that (also for / = k and a — ±) 

E«|P>±)(w±|P|<} = ^hlh*. 
Hence the claim follows by inserting an identity ((7j) and summing over k. □ 

3 Random dynamics of symplectic frames 

3.1 Symplectic frames and isotropic manifolds 

The space T v of symplectic p-frames, p = 1, . . . , L, is defined by 

T p = . . . ,u p ) \ ui e R 2L , (ui\u k ) = 6i >k , (ui\ J\u k ) — 0,l,k—l,...,p} . 

It is a manifold of dimension p(2L — p). One could also call T v an isotropic Stiefel manifold and 
J-l the Langrangian Stiefel manifold. 
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Proposition 1 The map C,:T L ^> 0(2L, R) n SP(2L, R) = XJ(L, C) de/med 6?/ 

CO) = (u,Ju) , u = (ui, . . . ,u L ) G Tl , 

is an isomorphism. 

Proof. This is immediate if one recalls 



0(2L,R) HSP(2L, 



a b 
—b a 



a,be M LxL (R) , a l a + b l b = 1 , a% = b l a 



Moreover, a + ib G U(L, C) gives the second isomorphism. □ 

A subspace £ C R 2i is called symplectic if (v\J\v') = for all v,v' G The isotropic 
manifold L p is by definition the set of all oriented symplectic p-dimensional planes in IR 2L . It 
is a manifold of dimension 2p(L — p) + \ p{p + 1). The maximal isotropic manifold L L is also 
called Lagrangian manifold. 

Next let A p C 2L , p = 1, . . . , L, be the vector spaces of the anti-symmetrized p-fold tensor 
products of C 2L . Decomposable (unentangled) vectors therein will be denoted by U\ A . . . A u p 
where ui G C 2L . A scalar product on A P C 2L is defined as usual by 

(ui A . . . A u p | u[ A . . . A u' p ) APC 2L = det p ((ui\u' k )i<i tk < p ) . 

As is well-known, oriented p-dimensional planes are isomorphic to the set of decomposable real 
unit vectors in A P C 2L . For the isotropic manifolds, this implies 

L p = {«! A . . . A v p G A P R 2L | ||ui A . . . A v p \\ = 1 , (vi\J\v k ) = , I, k = 1, . . . ,p} . 

Now each element u = (%,..., u p ) G T v defines a sequence of embedded, oriented symplectic 
planes u\ A . . . A u q G L q , q = 1, . . . ,p. Expressed in a different way, u G T p gives an element 
U\ A ... A Up G Lp as well as an unoriented, but ordered orthonormal basis therein. This is 
locally an isomorphism: 

Proposition 2 T v is a principal bundle over L p with fiber SO(p, R). 

The following elementary lemma about matrix elements of a Lagrangian projection w.r.t. 
eigenvectors of J (and hence also all rotations R constructed in Section ETTj) will be used later 
on. 

Lemma 2 Let IT = uu l be the projection on the Lagrangian plane associated with u G J-'l- If 
Vj G C 2L ; j = 1,2, are two normalized orthogonal eigenvectors of J, namely Jvj = tCjVj for 
signs o~j, then 

(1 + o-jCXfc) (vj\U\v k ) = 8 j)k . 
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Proof. Let LTi and LT2 be the projections on the first and second L components of C respec- 
tively. Hence LTi + n 2 = 1. Using the orthogonal (u, Ju) of Proposition [U one then has 

Sj t k = (vj\(u, Ju)(Hi + n 2 )(w, JuY\vk) = (uj-|im*|ufc) + (vj\ JuvfJ 1 ^) . 
As J* = J* = — J, the claim follows from the supposed properties of Vj. □ 



3.2 Action of a transfer matrix on a symplectic frame 

The group SP(2L,R) acts on the space Tl of symplectic L-frames. There is an obvious way 
to define such an action U : SP(2L, R) x Tl —> Fl (which will later on actually turn out to 
be relevant for the calculation of the Lyapunov exponents): given T gSP(2L,R) and u = 
(ui, . . . ,u L ) G Tl, the plane Tui A ... A Tu p is symplectic for any p < L; hence applying 
the Schmidt orthonormalization procedure to the sequence Tui, . . . ,Tul gives a new element 
Utu G Ti,. On a calculatory level, it will be convenient (and equivalent as one easily verifies) 
to define this action using wedge products: 

Ut Ul A ... A Up = —— —^rr , p = 1, . . . , L . 14 

\\A p Tui A ... A u p \\ 

More explicitly, this means that the pth vector of the new frame (Ut u) p satisfies for all v G C 2i : 

, ui, \\ ((A p_1 Tmi A ... A Up-i) A v\A p T u\ A ... A u p ) 

(v\(Ut u) v ) = —rn. — 77^ 1, I, . — rr— • (15) 

x IV JP/ ||Ap~ 1 Tm 1 A ... A Wp_i|| \\hPT Ui A ... A u p \\ V ; 

It is immediate from the definition that 



U ST = U S U T , S,T G SP(2L,M) . 

For fixed p, (fT4^1 defines a map Ut on L p . But the whole sequence p = 1, . . . , L, defines a map 
on Tl- Due to Proposition d this defines an action on the unitary group U(L, C). Therefore 
each Ut can be identified with an element of U(L, C) itself (explaining hence the notation with 
a letter U). 



3.3 Definition of random dynamics 

The random transfer matrices are symplectic so that they induce an action on the frames. Here 
their transformed normal form R(l — XP(n)) given by (J3J will be used to define a random dynam- 
ical system on Tl- The random orbits in Tl for some fixed initial condition (ui(0), . . . ,ul{0)) 
will be defined and denoted as follows: 



Ui(n) A ... A u L {n) = UR(i-xp( n ))Ui(n - 1) A ... A u L {n - 1) , n > 1 . (16) 

The free dynamics at A = is non-random and just given by the (elliptic-hyperbolic) rotation 
R, the analysis of which is straight-forward. The main objects of this work is to study the effect 
of weakly coupled randomness. 
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All the information needed for calculating averaged quantities like the Lyapunov exponent 
can actually be encoded in an invariant measure v on Tl- It can be defined by 

/ dv{u) f(u) = E / du(u) f{U R{l _ X p)u) ■ 

It is known |BLj that v is unique as soon as A 7^ and moreover continuous. In [JSS, Section 
4.7], perturbation theory for this measure (around the Lebesgue measure) was done in the case 
L = 1. The results below (in particular, Section FOJ) can be interpreted in a similar way. 



3.4 Channel weights 

It will be useful to introduce the probability (or weight) p p ,j(n) G [0, 1] for the pth frame vector 
u p (n) to be in the jth channel. More explicitly, 

ptM) = ( u p( n )Wt\u P (n)) , p pJ (n) = p+jin) + pp tj {n) . (17) 

The identity (J7J) and Lemma 121 (use 2vj eigenvectors of J to build a basis in 7TjC 2L ) imply 
respectively 

L c L 

^P P ,j{n) = !> J2p P j( n ) = v i ■ ( 18 ) 

j=0 p=l 

For an elliptic channel one has irj = tt^ so that 

Pp,j( n ) = 2 Pt,j( n ) ' N elliptic. (19) 

For the random dynamics of frames, let us introduce the mean presence probability of the pth 
frame vector in the jth channel: 

^ N-l N-l 

(Ppj)n = j^Yl E PpA n ) > (Pp,jP g ,k)N = E Ppd( n )P<iA n ) ■ ( 20 ) 

n=0 n=0 

Similarly, higher moments are defined. It follows from a standard ergodic argument that these 
quantities converge in the limit iV — ► 00 to some numbers denoted {p P j), {p P jP q ,k) and so on. 
However, this fact will not be used below. 



3.5 Separating hyperbolic from elliptic channels 

A bit of thought shows that the first 2L h + 1 frame vectors u 2 L h +i deterministically aline 

(up to an error) with the expanding hyperbolic basis vectors, that is w±, . . . , w^ h , w^_ L/ , . . . , w^ 
(recall that we supposed pi < 2, as in Fig. l(i)). Hence the remaining frame vectors have to be 
in the elliptic channels due to orthogonal and symplectic blocking. Hence let us call the frame 
vectors U2L h +2, ■ ■ ■ ,ul elliptic, while the first ones are called hyperbolic. The corresponding 
analysis, elementary but a bit tedious, is carried out below. It is not needed in case all channels 
are elliptic (as in Fig. l(ii)). 
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Proposition 3 For almost every initial condition (or almost every disorder configuration) , for 
n sufficiently large and for k = 1, . . . , Lh, 



|K|mW>| 2 = 1-0(X 2 ) , \(wt Aw+_ k \u 2k (n) Au 2k+1 (n))\ 2 = 1 - 0{\ 2 ) . 
Proof. Let us begin with the study of |(wj|wi(n)}| 2 . It follows from (fH^ and (fT^j) that 



e 2r!L |<w+| 


Mn))\ 2 


(ux(n)\R*R\ 


ui(n)) 



IMMn + 1)>|'= M — ^ + 0(A). (21) 
In order to analyse the denominator, let us insert (J7J): 



{ Ul (n)\R*R\ Ul (n)} = £ e^-^Kw+l^n))! 2 + e-^-^Kw^Mn))] 2 

= 1 + E (e 2 ^ - 1) |(w+| Ml (n))| 2 + (e^ - l)KwfMn)>| 



2 



Z=0 



where as before we identified indices L=0 and used that Ui(n) is normalized. The next aim is 
to prove an upper bound on this and therefore let us first note that e~ 2m — 1 < so that those 
terms can be discarded. For the remainder, starting from the smallest factor e 2m — 1 and going 
iteratively to the largest, one gets using each time \{w^\ui(n))\ 2 < 1 — X^=o I ( w t \ u i( n )) I 2 ' 



L h -l 



( Ul (n)\R*R\ Ul (n)) < e 2 ^ + ^ [e 2m - e 2 ^) |(w+Mr2))| 2 

(22) 



1=0 



< e 2 " 1 + (e 2 ^-e 2, ' 1 )|(w+| Ml (n))| 2 . 
Note that rji is the second largest hyperbolic angle satisfying 771 < Replacing in the above, 

,/ +. / u, 2 e 2 ^-^\(wj\ Ul (n)}\ 2 / \ \ 

w+K n + 1 2 > - 2( " 1V1/ +l > UI2 + ° A • 
1 + (e^ w — 1) \^^\u\\n))Y 

Thus, up to an error, \(w^\ui(n + 1))| 2 is larger than the image of x = |(w^|Mi(n))| 2 under the 
function f a (x) = 1+ ^L-^ X where a = e 2l - Vh ~' ni ' > 1. As n grows, this procedure is then iterated, 
giving rise to a discrete time dynamics through successive application of f a . The function f a 
has two fixed points in [0,1], an unstable one at and a stable one at 1. Either the initial 
condition is already away from the unstable fixed point or the random perturbation leads the 
discrete time dynamics to leave it (only with exponentially small probability one remains in its 
neighborhood, as an elementary argument shows). As it only takes a finite number of iterations 
to get within the neighborhood of the stable fixed point and the (random) perturbation is of 
order 0(A), one can conclude that, for n large enough, 
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|(w+K(n))| 2 = 1-0(A) . 

When this holds, however, the random perturbation cannot be linear in A anymore, because 
the 0(A)-term would not have definite sign and hence violate |(w^|«i(n))| 2 < 1. Indeed, it is 
elementary to verify also algebraically that the perturbative terms linear in A vanish in 
when one already knows (w^|wi(n)) = 1 — 0(A). Hence one can repeat the above argument in 
the neighborhood of the stable fixed point 1, but based on (}2*T]) with an error term 0(X 2 ). This 
implies the first claim. Moreover, due to normalization, | (w+|wi(n)) | 2 = 0(X 2 ) for all / ^ L as 
well as |(w^|?ii(n))| 2 = 0(A 2 ) for any I. 

The remaining estimates are proven by recurrence over k. After having exploited orthogonal 
and symplecting blocking, the basic argument is as before and therefore some calculatory details 
are suppressed and left to the reader. Hence let us suppose that | (w^ A w^_ ; |m 2 ; (n) A« 2 i+i (n) } | 2 = 
1 — 0(A 2 ) for all I < k. As {u2k{n)\u m {n)) = 5 2 k,m and (u2k{n)\J\u m (n)) = 0, one concludes 
that (u 2k (n)\w+) = 0(A) and (u 2k (n)\J\w+) = 0(A) for m = 1, . . . , k - 1, L - k + 1, . . . L. But 
for these hyperbolic channels J|w+) = |w~). As the same holds for u 2k +i(n), we can conclude 
that 



(w^ A wj\u 2k (n) A u 2k+ i(n)) = 0(A) , I, m = 1, . . . , k - 1, L - k + 1, . . . L . 
Using this, a short perturbative calculation starting from (|15jl shows 

1/ +a + I t ^i\a ( j_i\\\2 e 4 ^ |(w+ A wi_ k \ U2k (n) A ?i 2fc+1 (n))| 2 

\\ w k Aw l-k \u2k{n + l)Au2k+i{n + l))\ = -. r ~— ,,.„ ntnl 1—— -—- + 0(\). 

(u 2k (n) A u 2k +i{n)\A 2 R*R\u 2 k(n) A u 2k +i(n)) 



Generalizing the argument leading to (12211 . one can bound the denominator from above by 



e 



4vk+1 + (e 4% -e 4 ^)l(w^ Aw+_ k \u 2k (n) Au 2 k + i(n))\ 2 + 0(A 2 ) 



As r\k > r) k+ i, one can use the same function f a as above with a = e 4 ^ fc and complete the 
(two-stepped) argument as above. □ 



4 Lyapunov exponents 

4.1 Calculating Lyapunov exponents with symplectic frames 

In the definition (|2~Hj) of the Lyapunov exponents appears the operator norm. Instead, one may 
use symplectic planes as initial condition if an averaging over them is done. This is briefly 
discussed in this section. 

Important is the well-known fact |BL| that for any symplectic matrix T, its second quantized 
is most expansive on the isotropic subspaces, namely the norm of the second quantized A P T 
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can be calculated by ||A P T|| = sup ugip ||A p Tm|| = sup u&Fp ||A p Tw||. Furthermore, the Lyapunov 
exponents according to 1B1.1 A.III.3.4] the Lyapunov exponents are given by 



v 

1=1 



lim — E lo£ 



A' 



][[A P 2» u(0) 



n=l 



(23) 



APC 2 



where u(0) G T v is an arbitrary initial condition. One may average over u(0) w.r.t. to the 
invariant measure v on JF L and this immediately leads to 



7' 



du{u) E log ( || A P T u 



IApC 2L / ' 



where here the E is only an average over the single site transfer matrix T. Similar formulas can 
be found in |CL| Section IV. 6]. 



4.2 Basic perturbative formula: only symplectic channels 

As it is considerably more transparent, let us first perform the perturbative calculation of the 
Lyapunov exponents in the case where there are no hyperbolic channels (as in Fig. 1 (ii) ) . Hence 
we assume R to be orthogonal. 

Let us insert 1 = A p MA p M~ l in between each pair of transfer matrices in (}2*3*j) (the boundary 
terms do not change anything as can easily be argued as in |JSS| Section 4.1], for example). 
Then develop the product therein into a telescopic sum using the definition of the action (jnj) 
as well as the definition of the random dynamics of frames. This gives 



P j N-l 

= Yl Vlog(\\A p (M- 1 T(n + l)M)u 1 (n)A...Au p (n)\\ APC2L ) , 

1=1 n=0 

where E contains also an average over the initial condition it 1 (0) A. . . Au p (0). As R is orthogonal, 
so is A P R. Using A p (M~ 1 T(n)M) = A P R A p (l—XP(n)), one gets writing out the norm explicitly: 



^ = o J im m Z) Elog(det p (((l-AP(n + l))u,(7i)|(l-AP(n + l))u Jfc (7i))i<,, fc < p )) . 

1=1 n=0 

Now log detp = Tr p log, so that multiplying out gives: 

P N-l 

Y,H = o ,! im Tr E ETr p log(l p + (^H|(-A(P + P*) + A 2 |P| 2 )|n fc H) Wp ) , 

1=1 n=0 

where the argument n + lofP = P(n+l) was suppressed because they are all independent 
and identically distributed random variables over each of which can be averaged independently 
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in each summand. Finally let Tl p (n) be the projection in IR 2L onto the subspace spanned by 
ui(n),. . . , u p (n). Expanding the logarithm up to order C(A 3 ) and using that E Tr(H p (n)P) = 0, 
we obtain 



p A 2 1 / 1 \ 

= y Jjm - E Tr(n p H|P| 2 )--Tr(Pn»Pn») + 0(A 3 ) , (24) 

i=l n=0 ^ ' 

where the trace is now over IR 2L and P = P + P* is a real and self-adjoint matrix. Subtracting 
gives, up to C(A 3 ), 



7p = y Km - ^ E ( ( Mp (n)||P| 2 | Mp (n)) - ( Up (n)|Pn p (n)P|« p (n)) + -( Mp (n)|P| Mp H) 2 ) . 

n=0 ^ ' 

(25) 

As we shall show in the next section, the first and second contribution cancel exactly for the 
bottom Lyapunov exponent, while the third one can be calculated explicitly, namely it follows 
directly from Lemma El below (set = —1 therein so that there are no hyperbolic channels) 
that 

^ = i im ft E h 'M 2 - W (pljpl*)* + • ( 26 ) 

i,fc=0 

Because /i 2 > 1, Theoremnfollows immediately in the case where there are only elliptic channels. 



4.3 Oscillatory sums 

The aim of this section is to evaluate the terms appearing in the perturbative expansion (|25|) of 
the Lyapunov exponent. This will then give (|26[). As it does not take more effort at this point 
and will be needed below, we will however not suppose all channels to be elliptic, but only the 
exterior frame vector u p to be elliptic, i.e. p > 2L^ + 1. 

Lemma 3 Let p > 2L^ + 1 and suppose that the Main hypothesis holds. 
(i) 

, N-l L c 

-J2 V(u p {n)\\RP\ 2 \u p {n)) = - hi £ h] (p Pil ) N + 0(N-\\) . 

n=0 1=0 

(ii) 

, N-l j Lc 

-J] E( Mp H|P| Mp (n)) 2 = - ^ ^2(2-W OwW* + ©(AT*, A). 

n=0 j,k>L h 
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[Ill) 

, N-l L I L c \ L c 

E «n)\P\u q {n)f = 5 lE^? X>2<ft**>" + O^A). 

n=0 g=2L h +2 V Z>L h / fc=0 

Proof, (i) By inserting identities (|7j). 

(«p(n)| |i?P| 2 |«p(n)> = £ Mn)|(ir, + + Trf) |i?P| 2 (tt+ + tt* )M")> . (27) 

i,fc=0 

By Proposition El the sum may be restricted to Lh + 1 < l,k < L c at the cost of an error 0(A). 
Hence let us consider, for fixed elliptic channels /, k, and signs a, r, 



1 

J ( iV ) = a?E V MnM \RP\ 2 Trl\u p (n)) . 

n=0 

As from IfUjl . 

TIM")) = e™*<K(n-l)> + 0(A) , 
we get, because the boundary terms are of C(A^ _1 ), 

j(jV) = e <-°vi+TVk) j(]v) + 0(jV _1 ,A) . 

If now e l ^~ a " ni+Trik ' 1, this implies J(iV) = (^(iV™ 1 , A). By the main hypothesis this does not 
happen if I ^ k or a ^ fi. Therefore only the diagonal terms in (|27|) contribute to leading order 
so that 

. N-l L c N-l 

-J2 E(u p (n)\\RP\ 2 \u p (n)) = £ - £ £ E KWKI^I VM")> + O^A). 

n=0 l=L h +l n=0 cr=± 

Finally 7rf |u p (n)) = (| p P: i(n))^ \wf) for some complex unit vector u>f satisfying nf \wf) = \wf). 
Thus 

1 AT-l JV-1 L c 

^E e^wii^pw")) = atE E E E 2^ n )www + °( rl ' A )- 

n=0 n=0 Z=L h +l <r=± 

But the expectation value of the matrix element (over the random variable P only) is indepen- 
dent of wf and given by Lemma H^v). This directly leads to the first claim because the sum 
can again be extended to I = 0, . . . , L c by Proposition 01 

(ii) One has for p,q > 2Lh + 1 

(u p {n)\P\u q {n)) 2 = Yl Yl ( u p( n )K k PK l M n )) (u q {n)K T P^\u p {n)) . 

k,l,m,j=0 cr k ,cri,cr m ,aj=± 
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For the same reason as above, the sum can be restricted to elliptic channels k,l,m,j > Lh up 
to errors of order 0(A). From the 16 signs, Lemma ^iv) eliminates half, forcing = —&i and 
a m = —<Jj. To each of the finite number of remaining summands an oscillatory sum argument 
will now be applied. Set 

1 N ^ 

J ( N ) = A?E V(u p (n)\n-°Pn?\u q (n)) (u q {n)\^ P^\u p {n)) . 

n=0 

Proceeding as above shows 

J(N) = e 1(7(r,k+m) e iT( " m+ ^ ) J(N) + 0(N~\\) . 

Again invoking the main hypothesis, J(N) = 0(1) is therefore possible only if r = —a and 
{k,l} = {m,j}. Hence up to 0(N~ 1 , A), 

1 N ^ 

-J2 V(u p {n)\P\u q {n)) 2 

n=0 

^ N-l L c 

= nYI E E E ^)i^r^ri^w) KWkrpvrri^H) (28) 

71=0 k,l>L h <T=± 

N-l L c 

+ n E E E E {u P {n)K a P^i\<n)) (u q {n)KP^\u p {n)) . (29) 

n=0 k,l>L h ,k^l cr=± 

Normalizing the projections of the frame vectors and then applying Lemma ^iv) gives 

1 Lc ( 1 N ^ 

m = 2i E ti h i hvE eppMpM 

k,l>L h \ n=0 

The contribution (f2*9~|) can only be treated similarly if q = p. Supposing this, the second identity 
in Lemma H^iv) shows 

1 Lc ( 1 N ~ x 

(EHD = ]T hvE E feWfe( n ) 

k,l>L h ,k^l \ n=0 

The sum of the latter two contributions is given in (ii). 

(iii) One now has to sum (|28jl and ()29j) over g = 2Lh + 2, . . . , L, namely precisely the 
elliptic frame vectors. But because k and I only correspond to elliptic channels, the sum may 
be extended to q = 1, . . . , L because the weight of the hyperbolic frame vectors in the elliptic 
channels is of order 0(A) by Proposition El It now follows that the contribution of (J2T?|) vanishes. 
In order to show this, decompose 7r£ and 7rf therein using (fTTj) and note that the directions |w£) 
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and \wf) on which they project satisfy the hypothesis of Lemma El due to the identities (jlOj) . 
Therefore 

L 

^2 <KW> Mn)K = , k^l, (30) 

9=1 

implying the claim. The sum of (|28jl over g = 1 , . . . , L can easily be carried out using the 
identity (fTT^l : 

£ L c l c / N-l 

E® = ^ X>2 hvE 

<2=1 i>^h k>L h \ n=0 

Because w p is elliptic, the sum may carry over k = 0, . . . , L c because the error is 0(A 2 ). □ 
4.4 Sum of Lyapunov exponents near band center 

Let us again suppose in this section that there are only elliptic channels. Then it follows from 
(I24|) and Lemma El that 

L . 2 N-l L , v 

E 7i = y Mm - E E l P l 2 KM> " 2 («/H|Pn L (n)P|^(n)) + 0(X 3 ) 

1=1 n=0 i=l ^ ' 

iV— >oo o z— ^ 

i=l fc=0 

Using (|I8jl. one therefore gets: 

Theorem 2 Suppose that the Main hypothesis holds and all channels are elliptic. Then 

X> = ^ (/4) 2 + <?(A 3 ). (31) 
i=i 

Via the Thouless formula |CSj . this theorem also allows to deduce a 0(A 2 ) correction to 
the density of states. However, the term linear in A has to be calculated separately. Theorem 
12 allows to deduce an upper bound on jl, however, not a very tight one as we shall argue in 
Section El 

Corollary 1 Suppose that the Main hypothesis holds and all channels are elliptic. Then there 
exists a constant c such that 

1l < c A 2 log 2 (L) + 0(X 3 ) . 
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Proof. Because of the ordering of the Lyapunov exponents, it follows from Theorem El that 



2 



7, < | (H) 2 ■ 

But using sin(^) > a , one finds h\ v < ce log(L) for some energy dependent constant ce > 1- □ 

Finally, let us remark that it is also straight-forward to write out a perturbative formula for 
the top Lyapunov exponent: 



A 2 

7i = lim — 

N^oo 4 



Yl h ™ h ) (Pij)jv - 5^ ^( 2 - S j)k ) (pi,jPi, k )N 



2L 

i=o j,fe=o 



+ C(A 3 ). (32) 



This will be further analyzed in Section 



4.5 Weight of elliptic frame vectors in hyperbolic channels 

It follows from Proposition |3] and the arguments in its proof that the weight of a elliptic frame 
vector in the hyperbolic channels is of order 0(X 2 ), that is for N large enough 

(p Ptl ) N = 0{\ 2 ) , p = 2L h + 2,...,L , l = 0,...,L h . 

Actually, more detailed information about the leading order term as well as the redistribution 
of this weight on the contracting and expanding basis vectors will be needed below. As it turns 
out, the elliptic frame vector is randomly kicked into the hyperbolic channels and immediately 
forced back out; therefore it spends (to leading order) an equal amount of time in the expanding 
and contracting hyperbolic directions. 

Proposition 4 For p = 2Lh + 2, . . . , L and 1 = 0,..., Lh, one has 

(pti)N = (P^N + 0(N-\X 3 ) . 

Proof. Let us first calculate (wf\u p (n + 1)) in terms of (wf \u p (ri)) by using (jl5j) . In order to 
shorten the appearing expressions, let us drop the argument n in u p (n). The denominator in 
(|T3jl can be read off 



||A p J R(l-APKA...Au p || 2 = det((i?(l-AP)n fe | J R(l-AP)u m ) 1 < fcim < p ) = J] e 2 * + 0(X) 



k=0 



Therefore (wf\u p (n + 1)) is equal to 



Y[ e- 2r > k J (A P ^P(1 — AP) Ui A ... A «p_i A wf |A P P(1 — AP) % A ... A u p ) (1 + 0(A)). 
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The appearing scalar product in A P M. is given by 

/ (R(l - XP)u k \R(l - AP)u m ) 1 < fc/m<p _ 1 (R(l - XP)u k \R(l - XP)u p ) 1 < k < p ^ 1 \ 

det 

V (wf|i2(l-AP)Ol<m<p-l (Wf \R(1 ~ XP) Up ) J 

All the off-diagonal matrix elements are 0(A) due to the results of Section 13.51 except when 
a = +. In the latter case, the entries of the lower left corner are 0(1) for m = 21, 21 + 1 (again 
by Proposition |3J). Hence the contributions to the determinant up to order 0(A) are given 
by the product of the diagonal elements and (in the case a = +) by two transpositions. The 
diagonal elements of the upper left part are treated as above and cancel with the denominator. 
Therefore, 

(wfK(n + l)) = e ar » (wf|(l- XP)u p (n)) 

-5 a ,+e- 2r » J2 (™t\R(± ~ ^P>2i +J )(R(l - XP)u 2l+1 \R(l - XP)u p ) + 0(A 2 ) . 

3=0,1 

For the case wj~, one reads off 

(wf|up(n + l)) = e-"" (wf|(l- XP)u p (n)) + 0(A 2 ) , (33) 
while a bit of algebra invoking again Proposition El shows 

(w+\u p (n + l)) = e~* (w+\(l + XP*)u p (n)) + 0(A 2 ) , (34) 

Now set J a (N) = E — $2n.=o \( w i\ u p( n ))\ 2 - Because EP = 0, one gets by going back in 
history once 

1 N ^ 

J-(JV) = e-^J^(iV) + A 2 e-^E-^ \( W -\P\u p (n))\ 2 + 0(N'\ X 3 ) 

n=0 

1 1 ^ 

= A 2 ^ r - I E-^KHI^|w r )(w r |P| Mp (n)) + 0(iV- 1 ,A 3 ). 

n=0 

The appearing oscillatory sum can be treated by the same argument as in the proof of Lemma 
Ell). This gives 

J-(N) = A 2 J^^n^ X>* (Pp*)n + 0(N-\X 3 ) . 

k=0 

This argument can be repeated for J + (N) using (jBH) instead of (|HHJ). One finds J + (N) = 
J~(N) + 0{N~ 1 , A 3 ). Now the whole argument can be repeated for wj_ z . Finally summing the 
contributions of wf and wj_ ; allows to conclude the proof. □ 
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4.6 Perturbative formula for bottom Lyapunov exponent 

The aim is now to generalize the perturbative calculation of the Lyapunov exponents given in 
Section T4.2I to the case where there are both elliptic and hyperbolic channels. It is convenient 
to introduced the scaled frame vectors: 

ui(n) = e~ Vl ui(n) , fji = -(1 - 9 2 ^)V[l] • 

Note that, while the index on rji matches the channel index, the one on fji matches the frame 
vector: for a hyperbolic frame vector, fji is the expansion exponent in the direction into which 
u\ is alined by Proposition El but for an elliptic frame vector ui, one has fji = 0. Using the 
multilinearity of the determinant and then log det p = Tr p log, one finds 



p - N-l 

-Vi = ^ oiV ^ ^^p{^s{(ui(n)\ \R(l-XP)\ 2 \u k (n))i<i,k< P )) • 

1=1 ^ n=0 

The matrix elements of the leading order \R(1 — AP)| 2 = R*R + 0(A) now give a unit matrix 
l p . In fact, using the orthonormality property of the frame vectors and inserting (J7|), one finds 



({k{n)\H?R\u k (n))-5i, k = J] £ [e~^ e^'^ - l) («,(n)K|tt fc (n)) = 0(A). 

m=0 <7=± 

Moreover, if I = k this expression is 0(A 2 ). Indeed, for a hyperbolic frame vector ui, the 
summand m = I has a vanishing prefactor and all the others are 0(X 2 ) by Proposition El while 
for elliptic ui, all m corresponding to elliptic channels have vanishing prefactors and all the 
remaining hyperbolic m are 0(X 2 ) by Proposition 01 (actually, they were even calculated in 
Proposition 0J. Now around l p the logarithm can be expanded. In the expansion, the terms 
linear in P can be discarded because EP = 0. Hence Ym=i li ~ Vi is U P to 0(A 3 ) equal to 



7V-1 



n=0 L/=l 



l,k=l,tyk 



X 



2 P 



+ X 2 J2(Mn)\ \RP\ 2 \ui(n)) ~ — J2&i(n)\(R*RP + P*R*R)\Mn)Y 



i=i 



Lk=l 



When there are no hyperbolic channels, R is orthogonal and the formula reduces to ([25)1 . The 
bottom exponent can now be obtained by substraction. Let us suppose that fj^ — which 
means that E is in the spectrum of H L (0): 
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1L 



N-l 

lim — > E 



?1=0 



- {(u L {n)\R*R\u L {n))-l) 



1 i_1 



z=l 



A 2 



+ y (iti(ra)l |i?P| 2 |«z,(n)) 

y (u L (n)|(iri2P + P*i2*i2)| Ui (n)) 2 



X 2 



E (u L {n)\{R*RP + P*R*R)\ui{n))' 



1=1 



(35) 

(36) 

(37) 
(38) 

C(A 3 ) . (39) 



Now the terms (}35j) to (J39|) . each by definition containing the average ^^Lq 1 E but not 
the limit N — > oo, will be treated separately. Inserting (J7J) and using the normalization property 



m = \ E (e (1 ~^-l) (/>i)iv+(e-^^-l) (p^)* = E (cosh(2^)-l) H/) 

Z=0 Z=0 



iV 



the second step because the appearing averaged weights are equal by Proposition HJ Next, using 
orthogonality of u% and Ul, 



N-l L-l 



2 ^ E E E (EE (e* 1 "^ - l) KH|^|^(n)) 



n=0 J=l \k=0 o=± 

Again, the sum over k is actually restricted to the hyperbolic channels. But for hyperbolic a 
channel k, one has {uL{n)\7r%\ui(n)) = 0(X 2 ) unless a = + and I = 2k, 2k + 1 by Proposition El 
Hence 



AH 



\ E (^-l) 2 ^(pt fc )iv + 0(A 4 ) 



k=0 



which shows that to leading order and (pjj) compensate. The contribution (|3T|) was already 
calculated in Lemma Efi)- It will be compensated by ()39|) which is a bit more cumbersome to 
treat. Hence let us formulate it as a separate lemma. 

Lemma 4 Let p > 2Lh + 1 and suppose that the Main hypothesis holds. 



N-l 



^E E^KWKffflP + mil^n)) 2 = \hl f^hU Pp ,i)N + o(n-\\). 

n=0 1=1 1=0 
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Proof. Let us call the l.h.s. J(N). Because R*R\u k (n)) = e 2 ^\u k {n)) + 0{\), 

N-l L 

J ( iV ) = nYI E^( Mp (n)|(e^P + e^P*)|^H) (ui(n)\(e* P + P*)\u p (n)} + 0(A) . 

71=0 1 = 1 

Let us split J(N) = Jh(N) + J e (N) where Jh{N) contains the sum over indices 1 = 1,..., 2L h + l 
corresponding to hyperbolic frame vectors and J e (N) the remainder corresponding to elliptic 
frame vectors. In J h (N), one can replace up to 0(A) 

\ u i( n ))( u i( n )\ = *o > £ \ u 2k+j(n)}(u 2 k+j(n)\ = 7r+ , 

3=0,1 

by Proposition El Hence 

, N-l L h 

J ^ N ) = E^MnW^P + eKP^n+ieKP + e^P^lupin)) + 0(A), 

n=0 k=0 

and an oscillatory sum argument implies that Jh{N) is equal to 

EE E(va(e^P + e*P*)*+(e*P + e^P*)\vO + 0(N~\\) , 

m=0 cr=± k=0 

where w 1 ^ is some unit vector satisfying n^w^ = w a m . The expectation value of the last factor is 
by Lemma [U independent of w^. Now {pp^N = (Pp, m )v + 0(A) by Proposition 0] and (|19|). Of 
the four terms, the (P,P) and (P*,P*) pairs vanish after summing over a because of the sign 
in Lemma^ii). The remaining terms (P,P*) and (P*,P) are given by Lemma^iii) so that 

M*0 = 2 E h2 m(P a P ,m)» J E h l"k COsh ( 2 %) + 0{N~\X) • 
m=0 fc=0 

Finally for the / in the sum of J e (h), one has e^P + e _I *P* = P. Therefore, this is actually 
the term treated in Lemma HU^hi). Combining proves the lemma. □ 

Therefore, (|3*7j) and compensate to leading order just as do (|3*5|) to (|3l)]). The leading 
order contribution to 7l is thus solely given by (f5S]l. Noting that R* R\ul{ji)) = \ul{u)) +0(A), 
the contribution (|SH)l was already dealt with, in Lemma Efii) and we have proven: 

Theorem 3 Suppose that E e R «s m t/ie spectrum of Hl(0) and satisfies the Main hypothesis. 
Then 

n = J im £ h)hl(2-5 hk ) (p L , jPL , k )N + 0(A 3 ). (40) 

j,fc=0 
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The presented techniques also allow to write out formulas for the top Lyapunov exponent 
and the sum of the positive Lyapunov exponents. 

Proof of Theorem ^ (i) follows immediately from h\ > 1 and the fact that (pL,jPL,k) n is a 
probability distribution, (ii) For hyperbolic channels j, k = 0, . . . L^, the weights in ()40|) are 
O(X) by Proposition El so that they can be neglected. For the remaining elliptic channels j, one 

has h) > h\ c . But h 2 Le = l/sin(r/ Lc ) = (1 - ^Y l/2 = (e - <f )~ 1/2 > e" 1 / 2 . Replacing this 
concludes the proof. □ 



5 More insights on the channel weights 

This short section does not contain rigorous results. The aim is to get a better understanding of 
the averaged channel weights entering in the perturbative formulas above. Again, for simplicity, 
let us restrict ourselves to the situation where there are only elliptic channels. We first focus on 
the weights of the first frame vector U\. It follows from ()14|) that, up to 0(X 3 ), 



(ui|P*7r fe P|ui) - («i|7rjfc|wi)(wi| \P\ 2 \ui) 

+ (Mi|7Tfc|Wi)(Wi|-P|Mi) 2 - (Ui\{P*K k + ■K k P)\u 1 )(u 1 \P\u 1 ) , 

where the index n is left out on the r.h.s. and the expectation values is over P(n + 1) only. Now 
let us average E -h ^2 n =o anc ^ su PP ose that the limit exists. Then the l.h.s. vanishes. Hence 
the coefficient of A 2 on the r.h.s. has to vanish as well, up to O(X). Calculating it with an 
oscillatory sum argument as in Section 1431 shows that for all k — 0, ... , L c : 



Ep ljfe (n+1) -pi, k (n) = A 2 E 



L c , L 



= ttf ^2"khlh 2 (p hl ) - — ^2 u m h 2 m h 2 (p u p hk ) 



1L xr±i " 2L 

1=0 l,m=0 

+ 2l h 1 h m( 2 ~ S l,™) (Pl,lPl,mPl,k) ~ ^jr^2h 2 h 2 k (2-6 kjl ) (pi.jpi.fc). 

l,m=0 1=0 

These equations give relations between the averaged first, second and third moments of the 
weights py. Analogously, one can write out equations for (pijp^k) which then invoke up to the 
averaged sixth moments of the channel weights, and so on. This gives a hierarchy of equations 
for the channel weights. It results that the weights are independent of A and only depend on 
energy E (through the the frequencies if). 

It order to analyse these equations, let us close them already at first order by assuming 
factorization (pi,ipi,k) = (pi,i) (pi,k) and (pi,ipi, m pi,k) = (pi,i) (pi,m) (pi,k)- Furthermore we 
neglect the 5 m j and suppose v k = 2, both approximations which are 0(1/ L) w.r.t. the other 
terms. Now the sum over I factors and one obtains: 
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= hi - (pi,k) h ™ + (Phk) Yl h ™ (Phm) - h 2 k . 
m=0 m=0 

These equations have the unique solution (recall h\ = \j sin^)) 

1 + Z sm(%) 

where Z > is such that normalization Ylt^oiPhk) — 1 is assured. One easily verifies that 
Z ~ L. For small k (and large ones L — k as well) one has sin(^.) ~ k/L. Hence the weight 
on these channels is of order of unity, while it is of order 1/L on the others. But the channels 
with small and large k are precisely those near the band edges in Fig. l(ii) where the rotation 
frequency is small. Hence the weight of U\ is concentrated on the slowly rotating channels for 
which h\ = O(L) (they can be considered to be most similar to hyperbolic channels). Hence it 
is expected from (J3~2~|) that 7x = cA 2 + 0{\ 2 ) where c = 0(1) as L — > oo. Presumably, only the 
first few exponents are considerably larger than 7^. 

Due to symplectic blocking, the weight of U2 has to be centered on slightly faster rotating 
channels, eccetera. In conclusion, the weight of last frame vector ul is expected to be concen- 
trated on the channels which rotate the fastest and hence correspond to the band center. In 
these channels, h\ = 0(1) unless E is an internal band edge in which case h\ = 00. There- 
fore, away from these points the lower bound h 2 k >l which allowed to deduce Theorem^ from 
Theorem El is presumably not so bad because the weights of ul enter into formula pUj). In the 
case of mixed elliptic and hyperbolic channels, we expect the above argument to hold within 
the elliptic part of C 2i , namely for the weight vectors not alined to the hyperbolic channels by 
Proposition El 

Acknowledgments: This work profited from financial support of the SFB 288. While this 
paper was with the referees, we have done in collaboration with R. Romer extensive numerical 
studies of the perturbative formula (pfU|). It very well reproduces the energy dependence of the 
Lyapunov exponent (as calculated with the standard transfer matrix method), and this even 
for surprisingly large disorder strengths. 



References 

[AALR] E. Abrahams, P. W. Anderson, D. C. Licciardello, T. V. Ramakrishnan, Scaling Theory 
of Localization: Absence of Quantum Diffusion in Two Dimensions, Phys. Rev. Lett. 42, 
673-677 (1979). 

[APW] L. Arnold, G. Papanicolaou, V. Wihstutz, Asymptotic analysis of the Lyapunov expo- 
nent and rotation number of the random oscillator and applications, SIAM J. Appl. Math. 
46, 427-450 (1986). 

[Ben] C. W. J. Beenakker, Random-matrix theory of quantum transport, Rev. Mod. Phys. 69, 
731-808 (1997). 



27 



[BK] A. Bovier, A. Klein, Weak disorder expansion of the invariant measure for the one- 
dimensional Anderson model, J. Stat. Phys. 51, 501-517 (1988). 

[BL] P. Bougerol, J. Lacroix, Products of Random Matrices with Applications to Schrodinger 
Operators, (Birkhauser, Boston, 1985). 

[CK] M. Campanino, A. Klein, Anomalies in the one- dimensional Anderson model at weak 
disorder, Commun. Math. Phys. 130, 441-456 (1990). 

[CL] R. Carmona, J. Lacroix, Spectral Theory of Random Schrodinger Operators, (Birkhauser, 
Boston, 1990). 

[CS] W. Craig, B. Simon, Log Holder Continuity of the Integrated Density of States for Stochas- 
tic Jacobi Matrices, Commun. Math. Phys. 90, 207-218 (1983). 

[DK] H. von Dreifus; A. Klein, A new proof of localization in the Anderson tight binding model, 
Commun. Math. Phys. 124, 285-299 (1989). 

[FS] J. Frohlich, T. Spencer, Absence of diffusion in the Anderson tight binding model for large 
disorder or low energy, Commun. Math. Phys. 88, 151-184 (1983). 

[Fur] H. Furstenberg, Noncommuting random products, Trans. Amer. Math. Soc. 108, 377-428 
(1963). 

[GK] F. Germinet, A. Klein, Bootstrap Multiscale Analysis and Localization in Random Media, 
Commun. Math. Phys. 222, 415-448 (2001). 

[GM] Ya. I. Goldsheid, A. G. Margulis, Lyapunov exponents of a product of random matrices 
Russian Math. Surveys 44, 11-71 (1989). 

[GR] Y. Guivarc'h, A. Raugi, Proprietes de contraction d'un semi-groupe de matrices in- 
versibles. Coefficients de Liapunoff d'un produit de matrices alatoires independantes. Israel 
J. Math. 65, 165-196 (1989). 

[JSS] S. Jitomirskaya, H. Schulz-Baldes, G. Stolz, Derealization in random polymer models, 
Commun. Math. Phys. 233, 27-48 (2003). 

[KW] M. Kappus, F. Wegner, Anomaly in the band centre of the one-dimensional Anderson 
model, Z. Phys. B45, 15-21 (1981). 

[KLS] A. Klein, J. Lacroix, A. Speis, Localization for the Anderson Model on a Strip with 
Singular Potentials, J. Funct. Anal. 94, 135-155 (1990). 

[Klo] F. Klopp, Weak disorder localization and Lifshitz tails, Commun. Math. Phys. 232, 125- 
155 (2002). 

[PF] L. Pastur, A. Figotin, Spectra of Random and Almost-Periodic Operators, (Springer, 
Berlin, 1992). 



28 



[SSW] W. Schlag, C. Shubin, T. Wolff, Frequency concentration and localization lengths for 
the Anderson model at small disorders, Journal d'Analyse Math. 88, 173-220 (2002). 

[SB] H. Schulz-Baldes, Lifshitz tails for the ID Bernoulli- Anderson model, Markov Processes 
and Rel. Fields. 9, 795-803 (2004). 

[Tei] H. Teichert, Asymptotic analysis of the Lyapunov exponents of certain linear Hamiltonian 
systems with small noise, Math. Nachr. 172, 291-311 (1995). 

[Tho] D. J. Thouless, Maximum Metallic Resistance in Thin Wires, Phys. Rev. Lett. 39, 1167- 
1169 (1977). 



29 



